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Valleytronics in Tin (II) Sulfide 
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Tin (II) sulfide (SnS) is a layered mineral found in nature. In this paper, we study the two- 
dimensional form of this material using a combination of ah initio calculation and k • p theory. In 
particular, we address the valley properties and the optical selection rules of 2D SnS. Our study 
reveals SnS as an extraordinary material for valleytronics, where pairs of equivalent valleys are 
placed along two perpendicular axes, can be selected exclusively with linear polarized light, and can 
be separated using non-local electrical measurements. 

PACS numbers: 73.20.At, 73.43.Cd 


Introduction. The ever-expanding library of two- 
dimensional (2D) materials can be generally divided into 
several classes. One of them is composed of mono- 
elemental structures: the pioneering graphene and, as 
more recent additions, silicene and phosphorene. While 
graphene is the best known member of this group, phos¬ 
phorene has been experiencing a rapid growth in the sci¬ 
entific community due to a number of factors, including 
its mobility, high anisotropy, and gap tunability by strain 
and sample thickness. HHS] 

Another important class contains the transition metal 
dichalcogenides (TMDC’s), such as MoSe 2 , M 0 S 2 , WSe 2 , 
and WS 2 , among others. [6] These materials are semi¬ 
conductors with different gap sizes, a substantial spin- 
orbit coupling, and a valley degree of freedom [THS]- 
These properties make TMDC’s attractive for potential 
spintronic[T0] and valleytronic pTl applications. 

Recently, a new type of layered structures, referred to 
as group-IV monochalcogenides, has been brought to the 
attention of the physics community [13]. As the name 
suggests, these materials are composed of equal parts of 
group-IV metals and chalcogen atoms, and have a general 
chemical formula MC. In earlier works [T3UT7] . several 
of these materials have been investigated to determine 
their structural, electronic, and optical properties. These 
novel materials combine a number of features that make 
them very attractive for future study. Like TMDC’s, 
monochalcogenides also possess multiple valleys. How¬ 
ever, because of the rectangular unit cell, the valleys are 
situated on the axes of the Brillouin zone. Because of 
this, valley separation can be performed with linearly po¬ 
larized light instead of the circular polarization required 
for TMDC’s. In addition, perpendicular orientation of 
the valleys makes is possible to separate the valleys using 
the transverse non-linear conductivity. Finally, the pres¬ 
ence of heavy elements leads to a significant spin-orbit 
coupling. Thus, these materials are excellent canditates 
for optical, valleytronic, and spintronic applications. 

In this paper, we focus on one of the monochalco¬ 
genides: tin (II) sulfide (SnS). The rationale for our 
choice is the fact that SnS is found in nature in a mineral 
form, known as herzenbergite. To study the electronic 
properties of this material, we use ab initio calculations 


in conjunction with standard numerical and analytical 
methods. We develop a k • p model to describe the band 
structure at different points in the Brillouin zone and use 
these results to analyze the optical selection rules, as well 
as the nonlinear transverse conductivity due to the val¬ 
leys. Lastly, we discuss how those properties can be used 
to read and write the valley quantum number, laying a 
foundation for the use of SnS as a functional material for 
valleytronics. 

Electronic Structure. Monochalcogenides of group- 
IV elements are isoelectronic with phosphorus, and 
their structure is reminiscent of corrugated phosphorene. 
However, because of the two atomic species, the sym¬ 
metry of the crystal is lowered and bulk SnS belongs to 
the Pnma-D 2 ^ space group. Going from bulk to mono- 
layer, SnS loses the inversion symmetry, putting it in the 
Pmn2i space group. There are four symmetry transfor¬ 
mations in this group: the identity, a reflection, a re¬ 
flection with a translation, and a rotation with a trans¬ 
lation. The irreducible representations (irreps) for this 
space group, along with the operators belonging to each 
irrep are listed in Table [T| Following the standard con¬ 
vention, the axis of highest symmetry is labeled as z and 
it runs perpendicular to the corrugations. The ^-axis 
points along the corrugations; x axis is normal to the 
plane of the monolayer, see Fig. The direct products 
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TABLE 1. Irreducible representations and their transforma¬ 
tions under symmetry operations for the Pmn2i space group. 
Ri represents reflection across the plane normal to the i axis; 
Ci is the two-fold rotation around the i axis. The r before the 
operation means that the lattice is shifted by (0, ay/2, az/2) 
after the operation. 

of the Pmn2i irreps are listed in Table |H| 

The band structure of monochancogenides is distinct 
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TABLE IL Direct products of the irreducible representations 
for the Pmn2i space group. 



from both TMDC’s and phosphorene and contains two 
pairs of valleys, located on the ZT and FE axes. We 
obtain a few-band Hamiltonian at these valleys using the 
k • p formalism. 

The k • p Hamiltonian is constructed by finding the 
exact eigenstates and the corresponding eigenvalues at a 
point ko in the Brillouin zone, described by the Hamil¬ 
tonian Hq. One then treats the deviation from ko as a 
perturbation, given by 


Hi 


2 m 


+ — (p + ko) • k 

m 


(1) 


where Hq + Hi = iLfuii, m is the bare electron mass, p 
is the momentum operator, and k is measured from the 
position ko- 

We begin our analysis of the band structure by plot¬ 
ting the ab initio results along with the k • p approxi¬ 
mation in Fig. The first-principles calculations were 
performed using the Quantum ESPRESSO code. [ 18 ]. 
The exchange correlation energy was described by the 
generalized gradient approximation (GGA) using the 
PBEdl] functional. Interactions between valence elec¬ 
trons and ionic cores are described by Troullier-Martins 
pseudopotentials [ 20 |. The Kohn-Sham orbitals were ex¬ 
panded in a plane-wave basis with a cutoff energy of 
65 Ry, and for the charge density, a cutoff of 280 Ry was 
used. For the self-consistent calculation of the charge 
density, the Brillouin-zone (BZ) was sampled using a F- 
centered 1x10x10 grid following the scheme proposed 
by Monkhorst-Pack[2T]. For the calculation of the dipole 
matrix elements, a finer 1 x 40 x 40 grid was employed. 

Three independent k • p expansions are performed— 
one is around the F point, the other two are around the 
valley extrema. For this, we used 30 energy states with 
the corresponding matrix elements of Hi obtained from 
the density functional theory calculations. 

From Fig. it is clear that the k • p approximation 
is quite good at capturing the bands. Of course, using 
30 levels is impractical. Therefore, we construct smaller 
Hamiltonians around the centers of the valleys by ju¬ 
diciously choosing which bands can be dropped, while 
still describing the conduction and the valence bands cor¬ 
rectly. 

To obtain the k • p Hamiltonian, we calculate the ma¬ 
trix elements of Hi for the known eigenstates at ko, given 
by {i\Hi\j), where \i) and \j) are the eigenstates. It is 
clear from Eq. 0 that the terms proportional to and 


FIG. 1. (color online) k • p expansions around the F point 
and the valley extrema (color lines) plotted against the ab 
initio band structure (black lines). The crystal lattice of a 
monolayer SnS in two different projections and the Brillouin 
zone are also shown. 


ko • k are scalar and, thus, finite only if |i) = \j). On the 
other hand, the p-k = Pzkz-\-Pyky term can couple differ¬ 
ent bands. One can determine which matrix elements are 
finite based on the symmetries of the individual states. 

For a non-vanishing matrix element, the integral over 
all space of {i\Pn\j) has to be nonzero. This means that 
the term inside the integral must not be anti-symmetric 
for any of the transformations of the space group. In 
other words, the product must transform as the Ai irrep, 
see Table [H 

Let us first address the valley located along the ZT 
line. From the first principles, we determine that both 
the conduction and the valence bands belong to the Ai 
irrep. Starting with the matrix element for p^^ which is 
a part of the Ai irrep, we see that the direct product 
of the irreps for the bands and the operator is given by 
Ai^ Ai^ Ai = Ai, according to Table [H] leading to a 
finite matrix element for p ^. When we attempt to do the 
same for py {B2 irrep), we get Ai 0 B2 Z) Ai = B2 7^ ^1, 
resulting in a vanishing matrix element. In fact, both the 
valence and the conduction band couple only to B2 states 
in y direction {Ai 0 B2 Z) B2 = Ai). Thus, to ensure a 
correct dispersion in the y direction, we introduce more 
bands. It turns out that just two additional bands (the 
third valence and the third conduction bands, counting 
from those closest to the gap) are sufficient to capture 
the shape of the relevant valence and conduction bands, 
see Fig. |^a,b). The Hamiltonian at the ZT valley takes 
the following form 
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where £ zv is a matrix with the band energies 
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{Esv^ Ey, Ec, Esc) on the diagonal. The values of other 
matrix elements are provided in the Supplementary In¬ 
formation. 

Next, we move to the other pair of valleys, located on 
the YT line. Because the point of expansion is not on the 
high-symmetry 2 ; axis, the bands are no longer described 
by the four irreps listed above. Instead, the states fall 
into two categories A’ and A”, depending on their sym¬ 
metry across the yz plane. The coupling elements are 
finite only for the states in the same irrep, which is true 
for the valence and the conduction band. However, to 
obtain the correct dispersion, we need more bands. Just 
as before, the third valence and conduction bands are 
enough. The results for the YT valley are plotted in 
Fig. ic ,d). Here, the Hamiltonian is 
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The values of the interband coupling elements are also 
given in the Supplementary Information. 

Valley separation. Having determined the effective 
Hamiltonians for the two valleys, we are now in a po¬ 
sition to discuss the possibility of employing the valley 
degree of freedom. To do so, we look at the transitions 
between the tops of the valence and the bottoms of the 
conduction bands for the two valley pairs. The coupling 
strength with the light polarized in i direction is given 
by Ci = ^(c(k)||^|i;(k)). In other words, at k ^ 0, 
the coupling strength is proportional to the off-diagonal 
coefficient of the linear momentum term pointing in the 
direction of the light polarization. The transition rate 
is proportional to the coupling strength squared. As we 
established earlier, there is no coupling between the con¬ 
duction and the valence bands at the ZT valley in y direc¬ 
tion. This means that only light polarized in 2 ; direction 
can induce transitions between these two bands. At the 
TT valley, on the other hand, both y and 2 : polarized light 
can excite electrons. However, from the ab inition cal¬ 
culation, we see that the transition rate for y polarized 
light at this valley is approximately equal to the tran¬ 
sition rate for 2 : polarized light at the ZV. Also, both 
almost 40 times larger than the transition rate for the 2 ;- 
polarized light at the TT, see Supplementary Information 
for the matrix elements. This allows one to use linearly 
polarized light to selectively excite electrons from the va¬ 
lence to the conduction bands for individual valleys, see 
Fig. I^a). A similar finding was previously reported for 
TMDC’s, where the authors proposed using circularly 
polarized light for selective valley excitation. [zm 

Transverse Current. The final topic that we address is 



FIG. 2. Cross-sections of k • p results for a 4 x 4 Hamiltonian 
for (a,b) the ZV and (c,d) the YV valleys. The dashed lines 
are the k • p results, while the solid lines are obtained from ab 
initio calculations, (e) Conduction band profile as obtained 
from first principles calculations and (f) k-p conduction band 
profile at the minima, obtained from the k-p results for a 4 x 4 
Hamiltonian. 


the nonlinear response to the electric field in the context 
of multiple valleys. It can be shown that at zero temper¬ 
ature, the second order response to an applied electric 
field results in the following current 

C 

where r is the electron scattering time. This is a stan¬ 
dard result, but we derivation can be found in the Sup¬ 
plementary Information for the reader’s convenience. If 
the electric field points along one of the crystal axes (i), 
the transverse component of the integrand (pointing in j 
direction) becomes 

= (5) 

The current is nonzero only if the valley does not have a 
mirror symmetry in the j direction so the derivative with 
respect to qj does not give equal and opposite contribu¬ 
tion on the two sides of the valley. Using this reasoning, 
one can clearly see from the contour plot in Fig. [^that 
the transverse current is significant only for the pair of 
valleys perpendicular to the field. That is, if the field 
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FIG. 3. Valley selection and separation, (a) Valley selection 
for an external oscillating electric field. Depending on the 
polarization of the held, different valleys are excited, (b) Val¬ 
ley separation under a static electric held. Depending on the 
polarization of the held, diherent valleys how in the perpen¬ 
dicular direction. 

points in y direction, the transverse 2; current will only 
be observed for the ZT valleys. If, on the other hand, 
the held points in 2: direction, the transverse current will 
only arise in the VF valleys. 

Due to the valley structure, the transverse current in 
the individual valleys of every valley pair {ZT and VF) 


points in the opposite direction, so that the total val¬ 
ley current is zero. However, in a hnite system, different 
valley components accumulate on different sides of the 
sample (see Fig. [^b) for the current direction). This 
effect is similar to spin Hall effect, but instead of spin 
separation, here one has a separation of valleys. Thus, 
the valley population can be read through the applica¬ 
tion of an electrical current and measuring the transverse 
valley current using a valley hlter as described in Refill 
A similar device scheme can also be used to prepare a 
current with a dehned valley state. 

The analysis in this letter shows that the peculiar 
bandstructure of monolayer SnS and the rectangular 
shape of the underlying Brillouin Zone result in the pres¬ 
ence of the two pairs of valleys both in the conduction 
and valence band. Each pair of valleys can be optically 
pumped separately by excitation with linearly polarized 
light, thus writing the valley state. The valley state can 
be read using the reverse process, or alternatively by 
using the non-linear transverse valley conductivity aris¬ 
ing from the off-centricity and anisotropy of the valleys. 
This system thus allows for writing and reading without 
the need of circularly polarized light, or magnetic field, 
thus presenting clear advantages over other semiconduc¬ 
tor materials for spintronics and valleytronics. 
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